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Abstract 



Chiral symmetry serves as a guiding principle in low-energy hadron dynam- 
ics. An effective lagrangian, which explicitly breaks chiral symmetry via a 
small mass term, allows for a systematic method of calculating higher order 
corrections to tree-level results. This is applied to radiative muon capture on 
hydrogen, and to proton-proton 7r° production near threshold. Both processes 
' pose theoretical challenges in light of recent experimental results. 

(N ' 

^ ! I. INTRODUCTION 

m . 

lO . One of the more challenging research topics of modern nuclear physics is to relate the 

strong interaction dynamics of the quarks to the hadronic degrees of freedom of the nucleus. 
We know that the nucleons consists of at least three valence quarks, and we consider the 
quark dynamics to be determined by QCD. However, the predictions of QCD in hadronic 
low energy phenomena have proven too difficult to evaluate. For example, no one has been 
able to derive from QCD the strong nuclear forces between two nucleons. The analogous 
QED problem would be to calculate the van der Waals forces between two atoms, a task 
which has been accomplished successfully. Therefore, to investigate the low-energy strong 
hadronic interactions, the second best strategy would be to construct QCD inspired models 
(which all have their flaws), and/or to assume that the symmetry properties of QCD will 
manifest themselves also on the hadronic level. In this paper I will concentrate on hadronic 
manifestations and predictions of the almost perfect chiral symmetry of the QCD lagrangian. 
Provided the quarks are massless the QCD lagrangian is chiral symmetric. In this ideal case 
the chirality (or helicity) of a quark is conserved, meaning a right- and a left- handed quark 
will always retain their handedness. A quark mass term in the lagrangian will violate this 
symmetry and will allow left- and right- handed quarks to mix (or interact via spin-flip). 
Fortunately the u- and d- quark masses, m u pa 7 MeV and pa 15 MeV, are small compared 
to the QCD scale Aqcd ~ 200 MeV. It is therefore reasonable to consider the explicit 
breaking of chiral symmetry via the quark mass term in Cqcd as a small perturbation. (In 
the following we shall not consider the implied isospin violations which we know is small, 
but insist that m u = = m q .) From the Gell-Mann, Oakes and Renner (GOR) relation 

the pion mass squared: m\ oc m q . This means that the pions are massless in the chiral 
symmetric limit (when m q ^ 0). Below we shall argue that we can calculate the corrections 
to the chiral symmetry predictions of hadronic observables in a consistent, perturbative 
manner, via an approach called chiral perturbation theory (ChPT). 

After a short presentation of the essential points of ChPT, two examples of low energy 
hadronic physics will be presented where ChPT might give some valuable insight. Finally a 
brief overview of consequences of chiral symmetry considerations as applied to dense nuclear 
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matter and the question of possible kaon condensation in neutron stars. It was this last 
question and its possible consequences in supernova explosion by Brown and Bethe JTJ that 
triggered our investigations of ChPT. 



II. CHIRAL SYMMETRY AND THE EFFECTIVE LAGRANGIAN 

The seminal work of Gasser and Leutwyler has shown that ChPT is a very powerful 
and successful technique for gaining more insight of hadronic phenomenology at low energies. 
Its application includes the physics of the pseudoscalar mesons (ir,K,r)) and their interactions 
0, the pion-nucleon systems j| and to some extent the few-nucleon systems [|]||. 

The effective chiral lagrangian of ChPT that describes the low energy (E < A x ~ 1 
GeV) hadronic phenomena involves an SU(2) matrix U(x) which is non-linearly related to 



the pion field, e.g., U(x) = J 1 — [n(x)/ f n ] 2 + ir-7r(x)// 7r .[] In the meson sector, the sum of 

chiral-invariant monomials constructed from U(x) and its derivatives constitute the chiral- 
symmetric part of £ c h- The symmetry-breaking part of £ c h is given via a mass matrix Ai, 
the chiral transformation of which is dictated by that of the quark mass term in the QCD 
lagrangian. is written as an expansion in powers of and the mass matrix A4, both 

characterized by a generic pion momentum Q < A x , where Q represents either the pion 
momentum q or the pion mass m n . This means each term appearing in £ c h carry a factor 

(Q/A) u , where the chiral order index v is defined by v = d — 2, where d is the summed power 
of the derivative and the pion mass involved in this term. This suggests the possibility of 
describing low-energy hadronic phenomena in terms of £ c h that contains only a manageably 
limited number of terms of low chiral order. This is the basic idea of ChPT. 

The heavy fermion formalism (HFF) M allows us to easily extend ChPT to the meson- 
nucleon system. This HFF is used since Oq acting on the ordinary Dirac field ip describing 
the nucleon yields ~ M, a mass which is not small compared with A x . In HFF, ip is 
replaced by the heavy nucleon field N(x) and the accompanying "small field" n(x) through 
the "transformation" ip(x) = exp(—iMv ■ x) [N(x) + n(x)]. To obtain these new fields we 
use the projection operators V± = (l±^)/2 to project out the "heavy"-, N, and "small"-, 
n, fields from the Dirac spinor. By definition N(x) = V+exp(iMv ■ x)ip(x). Here the four- 
velocity is assumed to be static, i.e., ~ (1,0,0,0). A systematic elimination of the 
"small" field n(x) in favor of N(x) leads to an expansion in d^/M. Since M « 1 GeV ~ A x , 
an expansion in d^/M may be treated like the expansion in d fJi /A x . In HFF £ c h consists of 
chiral symmetric monomials constructed from U(x), N(x) and their derivatives in addition 
to the symmetry-breaking terms. Including the fermions the chiral order v in HFF is defined 
by v = d + n/2 — 2, where d is, as before, the summed power of the derivative and the pion 
mass, while n is the number of nucleon fields involved in a given term ||. As before, each 

term in C ch carry a factor (Q/A) u <^ 1. We use £ c h as given in M. 

In the most general form the effective lagrangian C c ^ involving pions and heavy nucleons 
in external weak- and electromagnetic-fields consistent with chiral symmetry, is written in 
increasing chiral order as: 



Here represents terms of chiral order u, and the explicit expressions for the C 
and C^Ni which only include terms of direct relevance for our calculations, are 



Another commonly used parameterization is U(x) = exp[ir-0(x)// 7r ]. 
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2 A 



tN l2M v ; 2M 2M 1 ' J+ 

+ 2cim£iViVTr([/ + W - 2) + (c 2 -^MiV(t;-u) 2 iV + c 3 Nu-uN 
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AM 

Above we have used standard notations 



DpU = d,U - i(V„ + AJU + iU^ - AJ; 

U = u 2 ; w M = iu^DfJJv) = i [u^d^u — ud^vf\ ; 



p R = v 4- A • F L — V — A 



Ft? = - d v F^ R - i[Fp' R , F^U /£ = + *F>t. (5) 

The covariant derivatives above include the external vector and axial vector fields, V M = V^ 2 ^ 
and A^ = AV^i respectively. The covariant spin operator of the heavy nucleon is defined 
by = 275^, J ft], which for our static velocity equals = (0, The coupling constants 
Cj,C2 and c 3 can be fixed from phenomenology Their values are related to the pion- 
nucleon cx-term, cj-kn (t) ~ (?/|m 9 ("U« + dd)[p) (where t = (jj — p) 2 ), the axial polarizability 
a a, and the isospin-even 7rA^ s-wave scattering length a + = |(ai/2 + 2a 3 / 2 ) ~ —0.008m" 1 . 
The four- Fermi non-derivative contact terms in Eq.(Q) were introduced by Weinberg || 
and further investigated in two- and three- nucleon systems by van Kolck et al. [Q. (In 
Eq.(g) 

rearrangement), only two of the four coupling constants C A are independent.) Although 
these terms are important in the chiral perturbative derivation of the short range nucleon- 
nucleon interactions not included in multi-pion exchanges @,[|, they do not play a major 
role in the following discussion. 

In addition to the chiral order index v defined for each term in £ c h , a chiral order index v 
is assigned for each irreducible Feynman diagram appearing in the chiral perturbation series 
M. Its definition for a multifermion system is 



the sum over A runs over the possible combinations of 7- and r- matrices: rf = 
= Sfj,, and Ty = S^t. However, because of the Fermi statistics (Fierz 



v 



A-E N -2C + 2L + Y.^ ( 6 ) 



where E^ is the number of nucleons in the Feynman diagram, L the number of loops, and C 
the number of disconnected parts of the diagram. The sum over i runs over all the vertices 
in the Feynman graph, and Vi is the chiral order of each vertex. One can show || that an 

irreducible diagram of chiral order v carries a factor (Q/A) u <C 1. However, the application 
of ChPT to nuclei involves some subtlety. As emphasized by Weinberg M, naive chiral 
counting fails for a nucleus, which is a loosely bound many-body system. This is because 
purely nucleonic intermediate states occurring in a nucleus can have very low excitation 
energies, which spoils the ordinary chiral counting. To avoid this difficulty, one must first 
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classify diagrams appearing in perturbation series into irreducible and reducible diagrams, 
according to whether or not a diagram is free from purely nucleonic intermediate states. 
Thus, in an irreducible diagram, every intermediate state contains at least one meson. The 
ChPT can be applied to the irreducible diagrams. The contribution of all the irreducible 
diagrams (up to a specified chiral order) is then to be used as an effective operator acting 
on the nucleonic Hilbert space. This second step allows us to incorporate the contributions 
of the reducible diagrams. We may refer to this two-step procedure as the nuclear chiral 
perturbation theory (nuclear ChPT). This method was first applied by Weinberg || and 
by van Kolck et al. |9| to the few nucleon system. Park, Min and Rho || applied the 
nuclear ChPT to meson exchange currents in nuclei, and others had success in describing 



the exchange currents for the electromagnetic and weak interactions (KJ[1T| 

In the literature the term "effective lagrangian" is often used to imply that the lagrangian 
is only meant for calculating tree diagrams. We must note, however, that the "modern" 
effective lagrangians have a different meaning. Not only can £ c h of Eq.(^) be used beyond 
tree approximation but, in fact, a consistent chiral counting demands inclusion of every loop 
diagram whose chiral order v is lower than or equal to the chiral order of interest. As will 
be discussed below, for a consistent ChPT treatment of the problem at hand, we therefore 
need to consider loop corrections and necessary lagrangian counterterms of chiral order v to 
cancel the loop-divergences. 

III. CHPT APPLICATIONS 

A. Radiative muon capture 

The radiative muon capture on the proton (RMC) pT +p — > n + v M + 7 has an extremely 
small branching ratio and to observe RMC is a great experimental challenge. Only recently 
has an experimental group at TRIUMF [|12|] finally succeeded in measuring the partial cap- 
ture rate Prmc(> 60MeV), corresponding to emission of a photon with > 60MeV. A 
main goal of measuring Trmc is to extract accurate information about the pseudoscalar 
form factor, fp, of the weak hadronic matrix element. The matrix element of the hadronic 
charged weak current h x = V x — A x between a proton and a neutron is 

A a\\ 



(n( Pf )\V x - A x \p(p t )) 



u 



(Pf) 



n ( 2\ A , /m(<? 2 ) \a , i 2\ A , fp(Q 2 ) A 

fv{q 7 + -7; a % + fA{q 7 75 + q 7s 



u(pi), (7) 



where q = pi — pf, and the absence of the second-class current is assumed. Ordinary muon 
capture on a proton (OMC), \i~ +p — > n + u^, gives only an approximate value for fp. The 
reason is that the momentum transfer in OMC, q 2 = — 0.88m^, is far away from the pion- 

pole position, q 2 = m 2 , where the contribution of fp{q 2 ) becomes most important. RMC 
provides a more sensitive probe of fp because in the RMC's three-body final state one can 
be closer to the pion-pole. Using the theoretical framework of Beder and Fearing the 
TRIUMF group extracted a value for fp about 50% larger than expected from PCAC. This 
surprising result should be contrasted with the fact that fp measured in OMC is consistent 
with the PCAC prediction within large experimental uncertainties. In this framework, as 
in many earlier works [jl4| , |l5fl , one invokes a minimal substitution to generate the RMC 
transition amplitude from the transition amplitude for OMC, the hadronic part of which is 

given by Eq. (|7]). Essentially, one writes fp(q 2 ) = fp/(q 2 — m 1)i where fp is a constant and 
then one replaces every q in Eq. ([?[) with q — eA (A is the electromagnetic field) except the 
q appearing in the q 2 dependence of /y, /a and fu- Due to the large discrepancy with the 

expected fp value, it seems reasonable to reexamine the reliability of this existing minimal 
substitution phenomenological approach ||13|| . 
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ChPT provides a systematic framework to describe the electromagnetic-, weak-, and 
strong-interaction vertices in a consistent manner, avoiding the phenomenological minimal- 
coupling substitution at the level of the transition amplitude. Furthermore, ChPT also 
satisfies the gauge-invariance and chiral-symmetry requirements in a transparent way. In 
the case of OMC, Bernard et al. [I6| and Fearing et al. [I7| used heavy-baryon ChPT to 



evaluate the value of fp with better accuracy but consistent with the value achieved in the 
PCAC approach. Muon capture is a favorable case for applying ChPT since momentum 
transfers involved here do not exceed m M , and is small compared to the chiral scale A ~ 
1 GeV, implying a reasonably rapid convergence of chiral expansion. 

A first calculation of the total capture rate r RMC and the spectrum of the emitted 
photons, dTnMc{k)/dk, in chiral perturbation expansion will be briefly discussed |]18||. We 



limit ourselves here to a next-to-leading chiral order (NLO) calculation and therefore we 
only consider the terms with v = and v — 1 in Eq. ([]]). To this chiral order we need 

only consider tree diagrams, and then C^ N only represents 1/M "nucleon recoil" corrections 

to the leading "static" part C^n- Furthermore, we do not consider explicit A degrees of 
freedom. The covariant derivatives of Eq.(|5|) include the external fields, V M and A^. These 
are determined by the W~ boson which couples to the leptonic and hadronic currents in 
the standard manner. In the actual calculation, taking the limit my/ — * oo, we make the 

1-2 

substitution — > (V~ — Aj l ) T ~ 2 lT , i.e. we treat V and A as static external sources. 
The only parameters appearing in the expressions for Eq.(^J) relevant for RMC [Eq.([|), first 
term of Eq.(|3|) and first and last line of Eq.(H)] are the pion decay constant, /„. = 93 MeV, 
the axial vector coupling, = 1.26, and the nucleon isoscalar and isovector anomalous 
magnetic moments, k s = —0.12 and k v = 3.71. Thus, to the chiral order of our interest £ c h 
is well determined. 

We consider all possible Feynman diagrams Figs.|l]-^| up to chiral order v = 1 which 
contribute to the process /i" + p — > n + u + ^. The zigzag lines in these diagrams repre- 
sent the W~ boson. For static W~ bosons the diagrams in Figs.[l]-|6] reduce to those that 
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would result from the simple current- current interaction of the V — A form. The reason for 
explicitly retaining the W~ boson lines is to clearly separate the different photon vertices 
(see e.g. Fig||). The leptonic vertices in these Feynman diagrams are of course well known. 
The hadronic vertices are obtained by expanding the ChPT lagrangian [Eqs. (jl]), @, © 
and the relevant parts of (f|)] in terms of the elementary fields N, tt, V and A and their 
derivatives. The evaluation of the transition amplitudes corresponding to these Feynman 
diagrams is straightforward. We denote by Mj (i = 1 . . . 6) the invariant transition ampli- 
tudes corresponding to Fig. ([!])- (|6|), respectively. They are calculated in the Coulomb gauge, 
e(A) = 0, and are given by: 



i.e. v 



M 1 = e /3 (A) 



- / \ /-i \lr ft + m M / i\ 

U V (S)~{ T {1 - 75 ) IfjU^S ) 



2{k- fi) 

M t = [u v {s)lr{l-lsW{s')] \Nl(a)hl(X)N p (a r 



Ni{a)K[N p {o>) 



(8) 



2,3,4,5,6 



where hj; i 



1, ■ ■ • , 6 are the hadronic operators given in Ref. flSfl . In Eqs.(|) and @ fi, 

v, p = (E p ,p), n = (E n ,n) and k = (cu^, k) are the four-momenta of the muon, neutrino, 
proton, neutron and photon, respectively. The z-components of the spins of the muon, 
neutrino, proton and neutron are denoted by s, s', a 1 and a, respectively, while e(A) stands 
for the photon polarization vector. 

We consider only the case of RMC from the /i-p atom with statistical spin distributions, 
leaving out the hyperfine-state decomposition and the treatment of RMC from the pfip 
molecule. With our kinematical approximations: both the muon and the proton are at rest 
and neglecting the recoil neutron kinetic energy, the spin-averaged total capture rate is 



RMC 




|$(0)p 



(2tt)< 



d 3 



n 



x5^\n + v + k — p 



(2tt) 3 
A*) E 

(tct'ss'A 



d 3 



V 



(2tt) 3 
|M| 2 



d 3 k 



(2tt) 3 2oo k 



(9) 



m^. The numerical values for the total capture rate r RMC are then |L8 



Trmc = 0.061 s 1 (excluding the pion-pole diagrams' con- 



where the sum is over all spin and polarization orientations, M = X^=i Mi, and the \ip 
atomic wavefunction at the origin is $(0). Within the kinematical approximation stated 

earlier, 5^(n + v + k — p — fi) w 6(\k\ + \ v\ — m^S^ (ft + V + k) and the maximal 7 energy 

IS (^fc)max 

1. To lowest chiral order v - ' ■ 1 n 
tributions Trmc = 0.043 

2. Including the v = 1 (0(1/M) recoil terms) contributions give Trmc 
pion-pole = 0.053 

Our total capture rate Trmc = 0.075s -1 is close to the value given in JTJJ], Trmc = 0.069s -1 , 
and practically identical to Trmc = 0.076s _1 reported in |15|]. Our 0(1/ M) recoil correc- 
tions account for about 20% of the leading order 0((1/M) ) contribution, which indicates a 
reasonable convergence of the chiral expansion. As one can see about 30% of the total value 
of Trmc comes from the pion-pole exchange diagrams in our calculation. 



0.075 s 



-1 



no 



A direct comparison of our calculation with the experimental data | |12| is premature 
because we have not considered captures from the singlet and triplet hyperfine states sep- 
arately, or capture from the p\xp molecular state. This also means that at this stage we 



cannot directly address the u fp problem" that arose from the TRIUMF data |12j| . However, 
as already stated ChPT gives a unique prediction on fp, are consistent with theGoldberger- 
Treiman value fp = Q.QgA fli6Ul7| . Meanwhile, for the spin-averaged /zp-atomic RMC, our 
ChPT calculation gives a photon spectrum that is harder than that of |H| for the PCAC 
value of fp. Of course, a more quantitative statement can be made only after a more 
detailed ChPT calculation becomes available in which the hyperfine states are separated 
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and the p/xp-molecular absorption is evaluated. We must also emphasize that the present 
calculation includes only up to the next-to-leading chiral order (NLO) contributions. A 

next-to-next-to-leading order (NNLO) calculation includes the v = 2 chiral lagrangian, 

and C^jf, and also loop corrections arising from and has been completed and these 
next order terms are found to give small corrections to the tree diagrams (The fi- 

nite contributions from the loop diagrams would give momentum-dependent vertices, which 
would correspond to the form factors in the language of the phenomenological approach 
13l , |l5ll .) We note that the approach of Bernard et al. 



which we have used here, does not 
contain the explicit A degree of freedom in contrast to p0| . For a complete calculation the 
A degrees of freedom should be included in the next chiral order as done by, e.g., Fearing 
et al. P7|. 



B. Pion production at threshold 

The second application of ChPT is concerned with recent high-precision measurements 
of the total cross sections near threshold for the reaction p + p — > p + p + ir° by Meyer et al. 
pl| . Traditionally the pion production reactions are described by the single nucleon process 
(the Born term), Fig.7(a), and the s-wave pion rescattering process, Fig. 7(b). 




Pi p[ 



— p p — 



— p p. 




Pi p[ 



The Born f errri is assumed to be given by the pseu§ovector interaction Hamiltonian 

H ° = V ( T-,r ) - 2lf { °" V ' T '^ } ) (10) 

where Qa is the axial coupling constant, and f n = 93 MeV is the pion decay constant. The 
first term represents p-wave pion-nucleon coupling, while the second term accounts for the 
nucleon recoil effect and s-wave pion coupling. The s-wave rescattering vertex in Fig. 7(b) is 



commonly calculated using the phenomenological Hamiltonian |2~2"| 



T~Ci = All ——IplT ■ TT1p -\- AlT——-lfjT- 7T X 7T1p (11) 

The two coupling constants Ai and A2 determined from s-wave pion nucleon scattering 
lengths have the values Ai ~ 0.005 and A 2 ~ 0.05. Thus, Ai <ic A 2 as expected from 
current algebra. The early calculations f22|-p4|j which are based on these phenomenological 
vertices, underestimate the measured 7r° production cross sections by a factor of ~5. Since 
the second term in Eq.([lT]) is suppressed compared to the first term by ~ m n /M, the 
importance of the rescattering term is enhanced. However, A 2 is much larger than Ai, and the 
isospin structure of the A 2 term is such that it cannot contribute to the ir° production from 
two protons at the rescattering vertex in Fig. 7(b). Thus, the use of the phenomenological 



S 



Hamiltonians, Eqs.(fL0|) and (|TT|), to calculate the Born term and the rescattering terms, gives 
significantly suppressed cross sections for the pp — > ppn reaction near threshold compared to 
7r^ production because only the small correction terms contributes. Therefore, the calculated 
cross sections can be highly sensitive to any deviations from this conventional treatment. 

It is convenient for our discussion to introduce the typical threshold kinematics for this re- 
action. Consider Fig. 7(b) in the center-of-mass system with the initial and final interactions 
turned off (since it will modify what follows). At threshold, (q , q) = (m^, 0), p' 1Q = p' 20 = M, 
and p[ = p' 2 = 0, so that any exchanged particle must have ko = m n /2 = 70 MeV and 



= ym n M + (m 7T /2) 2 ~ 370 MeV/c, which implies k 2 = —m^M. Thus the rescatter- 
ing process probes internucleon distances ~ 0.5 fm. The reaction is therefore sensitive to 
exchange of heavy mesons which are important in phenomenological N-N potentials. Lee 
and Riska |25 showed that heavy- meson exchanges (scalar and vector exchange) could be 



capable of enhancing the cross section significantly. Meanwhile, Hernandez and Oset |26 



and Hanhart et al. p7j showed that the off-shell dependence of the nN s-wave isoscalar 
amplitude featuring in the rescattering process, k 2 = —m^M ^ m 2 , could also appreciably 
enhance the rescattering amplitude. Given these developments we consider it important to 
examine the significance of these phenomenological lagrangians in ChPT. 

Systematic studies based on ChPT would be valuable to sharpen our conclusions on 
whether or not the heavy-meson exchange contributions are needed to explain the observed 
cross section f or p p — > ppir . Three very similar ChPT investigations have recently been 
completed p8|-|30[|. Here we describe briefly the work of the USC group. To generate the 
one- and two-body diagrams of Figs. 7(a) and 7(b) we minimally need terms with v = 

and 1. Eq.(|I]) leads to the pion-nucleon interaction Hamiltonian Hi nt = TQnt + ?4nt' w ith 
Wi2 = If N[^^(r-7T)]N + ^-NT-7zxnN (12a) 



«2 = —^N{<r-V,T-n}N + l[2 Cl m^ 2 -(c 2 --^)7r 2 -c 3 (97r) 2 ]iViV (12b) 

Here Ti^ t represents the term of chiral order i>. We now compare Timt resulting from ChPT, 
Eq . (|T2"D with the phenomenological effective Hamiltonian TLq +Hi, Eqs.fllOD and (]TT|) . We 

note that the first term in H.^ and the first term in TiS 1 ' exactly reproduces Ho- Thus the 
so-called "Galilean-invariance" term naturally arises as a 1/M correction term in HFF. As 

for the ttttNN vertex, we can associate the second term in Tif^t to the A2 term in Tii, and 

second term in Ti^t to the Ai term in Hi. Consequently, 

4nX-/m n = f [2 Cl - (c 2 - - c 3 g] ee n(k,q) (13) 

with q = (u q , q) and k = (u^, k). Now, for on-shell low energy pion-nucleon scattering, i.e., 
fc~g~(m vr ,0), we equate 

47rA^/m 7r = /t = K(k = (m^,0),q=(m^,0)), (14) 

where 

ml ( g\\ ( m n \ 3g 2 A ml 

K0 = 7 2 [ 2Cl - C ^ C3+ 8M) = - 27V { 1 + ^) a + 12^T/f (15) 

The above cited empirical value for a + leads to k q = (0.87 ± 0.20) GeV -1 . If we keep only 
the lowest chiral order, the first term in Eq . (|i~5|) , A* is 
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47rA^ 



= -2tt ( 1 + 



m N 



= (0.43 ± 0.20) GeV" 1 , 



(16) 



or A* = 0.005 ± 0.002, identical to the "standard value" of Ai used in the literature, 
Eg. (|TTD . We note that the ChPT value for k = 0.87 GeV -1 , Eq.(|T3j), include next chiral 
order terms oc g A , and is twice the value of the first term, Eq.(|16|). In our comparison 
between the traditional and the ChPT approaches below, we shall use Eg. (|T3|) . Obviously, 
since n(k, q) depends on the four-momenta q and k, we cannot identify k with the constant 
X\. To illustrate how the q and k dependencies in n(k,q) affect the rescattering amplitude 
of Fig. 7(b), we consider the typical threshold kinematics, g~(m 7I 
and denote this n(k, q) value by Kth- 



0) and fc~(|m T , ^m^M), 



Kth 



J 7T 



2 1 8m N 2 



■1.5GeV" 



(17) 



Thus the s-wave pion-nucleon interaction is not only much stronger than the on-shell cases, 
of Egs.(|l5D and (0), but the sign of the off-shell coupling strength is opposite to the on-shell 



cases. The first feature is gualitatively in line with the observation of Refs. [p6|,27l that 



the rescattering term should be larger than previously considered. However, the sign of the 



typical off-shell coupling, n t h, is opposite to the one used in Refs. |26],|27[. This flip of the sign 



drastically changes the pattern of interplay between the Born ana rescattering terms. The 
sign change in ChPT arises from the ttN rescattering vertex where the irN scattering occurs 
at an energy below the nN threshold. If we force the exchanged pion on the mass-shell and 
only consider the exchanged pion three- momentum \k\ to be off-shell, we essentially recover 
the results of Hanhart et al. [27] . 

The two-nucleon transition matrix element T for the pp — > ppir process is T = 
($/|T|$j), where (|$/)) is the initial (final) two-nucleon state distorted by the initial- 
state (final-state) interactions. Here T represent the contributions of all irreducible diagrams 
(up to a specified chiral order v) for the pp — > ppir process, and we use T as an effective 
transition operator in the Hilbert space of nuclear wavefunctions. 



T = T (_1) + T {1) 



>T-Born I >T-Res 
1 -\ 2 +l 



(18a) 



✓T-Born 
1 -\ 



9a 



Am N f n 



U), 



i=l,2 



(18b) 



rr Res 
2 +l 



9a ^ (l , a-fkiif 

J tv i=1 2 K i m n + %e 



(18c) 



The operator T^ v > represents the contribution from Feynman diagrams of chiral order v. 
The initial and final three- momenta of the i-th proton are Pi and p\, hi = pi — p\\ and 
n(ki,q) is as defined in Eq.(^). To chiral order v = 1 we have additional contributions 
due to the loop corrections. These loop corrections to the Born diagram, Fig. 7(a) introduce 
an effective form factor at the vertex defined by C c h- The loop corrections do not change 
drastically our results. Since we cannot reproduce the measured cross section, we ignore 
these corrections in this paper. 

A formally "consistent" treatment of T would consist in using for |$,) and |$/) two- 
nucleon wave functions generated by irreducible diagrams of order up to v — 1. A problem in 
this "consistent" ChPT approach is that the intermediate two-nucleon propagators in Fig. 7 
can be significantly off-mass-shell. Another more practical problem is that, if we include 
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the initial- and final- two-nucleon (N-N) interactions in diagrams up to chiral order v = 1, 
these N-N interactions are not realistic enough to reproduce the known N-N observables. A 
pragmatic remedy for these problems is to use a phenomenological N-N potential to generate 
the distorted N-N wavefunctions. Park, Min and Rho |TT[ used this hybrid approach to 
study the exchange-current in the n+p — > r y + d reaction ana at least, for the low-momentum 



transfer process studied in Ref. |[L1|| , the hybrid method is known to work extremely well. 

Following standard practice in nuclear physics the T-matrix was evaluated using r-space 
wave functions after Fourier transforming T B ° rn and T5 es , Eq . ( 18 ) , into r-space. To simplify 



the Fourier transformation the fixed kinematics approximation, K,(k,q) = Kth, Eq. (fL7D , was 
used |28| , 29[| . The results of these two groups indicate that, for the various nuclear distortion 



potentialsconsidered, the calculated cross section is much too small to reproduce the experi- 
mental cross section. If we define the discrepancy ratio R by R = a^t '/ 1 ^tot \ with o^f taken 
from Ref. jnj], then R = 80 (R = 210) for the Hamada- Johnston (Reid soft-core) potential, 
and R happens to be almost constant for the whole range of Ef < 23 MeV for which o^t is 
known [p9[| . Thus, the use of Kth, Eq . (|l7l) , results in a significant cancellation between the 
Born- and the rescattering terms. Thisdestructive interference is in sharp contrast to the 
constructive interference of Refs. [22|,|25|-p7[] . However, using ChPT in a systematic fashion 
we have shown that the contribution of the pion rescattering term can be much larger than 
obtained in the traditional phenomenological calculations. 

We also learn that the fixed kinematics approximation (which is commonly used in the 
literature) should be avoided. There are at least two reasons why this is not a good ap- 
proximation for this reaction: (i) The initial- and final-state interactions play an essential 
role in the near-threshold pp — > ppn reaction; (ii) The theoretical cross section within the 
framework of the Born plus rescattering terms is likely to depend on the delicate cancellation 
between these two terms. In a just completed momentum space calculation we avoid the 
n(k,q) = K t h and the fixed kinematics approximation |j30"| . We use Eq.fllBD directly where 



K(k,q) is given by Eq. (|T3"D . Our results show that the magnitude of T5 es increases by a 
factor ~ 3 which means the irN rescattering term, Fig. 7(b) dominates |3(J. Further, we 
confirm the sign of the rescattering amplitude found in the first two ChPT calculations 
f2l||2"9"|| , and we find that the c\ term in K(k,q) dominates since the irN rescattering terms 
C2 and C3, Eq.(|T3"|) average to a tiny value in the distorted wave integrals due to their energy 
dependences. The constant c\ is given by the nN sigma term, E„-Ar(0), and the uncertainty 
in the magnitude of the rescattering amplitude, T+f 3 , is given directly by the uncertainty 
in the value of X^^O). Due to the destructive interference with the Born term, this un- 
certainty is enhanced in the resulting cross section. However, our calculated cross section 
is still too small. Presently we |3]J are evaluating the next chiral order two-pion exchange 
diagrams which simulate partially the cx-meson exchange used in Ref. [p5[] . To achieve a large 
cross section these terms should be important and this call into question the convergence of 
ChPT. Phenomeno logically one could also include the repulsive u;-exchange when we com- 
pare with experiments. In the effective C c h of Eq.(|3[) this is included in the four-nucleon 
terms, where we would calculate the constant Ca as given by us exchange. For charged pion 
production from pp, the higher chiral order terms considered here are smaller corrections to 
the dominant Born- and rescattering terms. For 7r° production, however, these next order 
chiral correction terms give the dominant amplitudes. 

Finally, we should mention that when ChPT is extended to SUp(3) and we assume the 
KN sigma term Exjv(O) ~ c\ to dominate, the effective Kaon mass in nuclear matter could 
become very small and we could have a Kaon-condensate in dense nuclear matter. However, 
if we use current algebra, "weak" PCAC and standard nuclear physics treatment of meson 
self-energy in matter, we do not find a Kaon condensate in matter. This is still an open 
question and the existence of a Kaon condensate in dense nuclear matter opens exciting 
possibilities in astrophysics, e.g. Ref. M. 
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